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Abstract 

We study a new iP —>■ tttt decay channel with gluons in intermediate state which is 
normally neglected within the factorization framework. Both short-distance and long¬ 
distance parts of the amplitude are calculated. The chiral Lagrangian approach is used 
for obtaining the long-distance contribution. The nonperturbative contribution gives an 
additional enhancement to K —>■ tttt decay amplitude with AJ = 1/2. A sizable violation 
of factorization in the p'^ order of chiral perturbation theory is found. 
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The origin of considerable enhancement of AJ = 1/2 parts of non-leptonic kaon decay 
amplitudes remains one of subtle points within the standard model (SM) [1]. Though ex¬ 
plained qualitatively by strong interaction effects it still escapes the reliable quantitative 
description. Numerous attempts have been recently made to achieve a sufficient accuracy 
for non-leptonic kaon decays in SM. The efforts were directed to improving the pertur¬ 
bative QCD analysis [2] and to accounting for long-distance effects using more advanced 
models of strong interactions at low energy, for example, the chiral Lagrangians [3], l/N^ 
expansion of QCD [4,5] or lattice simulations [6]. The present results however do not £t 
experimental data. 

An effective AA = 1 Hamiltonian reads [7-9] 

HAs=i = ^VudV*, Y. [ziM + ^yiif^)]Qi ( 1 ) 

where Gp is a Fermi constant, V is the quark flavor mixing matrix, r = —VtdV^g/VudV*^, 
Zi^jj) and are the coefficients of Wilson expansion, {Qi} is a full basis of dimension 
six local operators containing light quark helds (m, d, s) only. 

The renormalization group improved perturbation theory does not account for strong 
interaction of soft light quarks and gluons. The information about this interaction is 
entirely contained in the hadronic matrix elements of local quark operators. Factorization 
procedure for evaluation of these matrix elements [10], i.e. the procedure of replacing the 
four-quark operators by a product of two non-interacting quark currents accounts only 
for the ’’factorizable” part of this interaction [11]. But there are also ” non-factorizable” 
contributions, for example, ones corresponding to annihilation of a quark pair from the 
four-quark operator into soft gluons, which are omitted within the factorization procedure. 
The calculation of these contributions and the generalization of matrix element estimates 
beyond the factorization framework can be systematically done within QCD sum rules 
technique combined with the chiral Lagrangian approach. This possibility is related to 
studying a new K tttt decay channel generated by annihilation of a quark pair from the 
four-quark operator into gluons with the subsequent formation of the final pion pair by 
the soft gluon cloud, i.e. the decay channel with gluons playing the role of an intermediate 
state. Being non-factorizable, this decay mode does not appear as a correction to some 
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leading order contribution and can be studied by its own. This feature makes obtained 
results more accurate. 

In the present paper we study a new K tttt decay channel with the simplest scalar 
colorless gluon conhguration forming an intermediate state. We calculate both short- 
distance (perturbative) and long-distance (nonperturbative) parts of the corresponding 
amplitude. To obtain the long-distance contribution, the chiral effective Lagrangian is 
used as a low energy model of strong interactions and the relevant factor of proportionality 
between quark and mesonic operators is derived via QCD sum rules. 

Non-factorizable contributions reveal themselves in two different ways: hrst, they 
appear as corrections to the couplings characterizing the ’’factorizable” weak chiral La¬ 
grangian in and higher orders, second, some new non-factorizable terms emerge. 

The latter is the case for K tttt decay mode with gluons forming an intermediate state. 

Before turning to the long-distance effects of meson-gluon transitions it is useful to 
consider the similar phenomenon arising already in perturbative QCD as a leading correc¬ 
tion in the inverse mass of charmed quark. The effective low energy tree-level Hamiltonian 
for AS = 1 transitions before decoupling of the c-quark reads 

= ^K.,v;;(Q 5 - (i - r)Qg (2) 

where Ql = 4(si:,7^gL)(gi7^(ii;,), q = u,c and qL(R) stands for the left(right) handed 
quark. Performing the OPE to the hrst order terms in a* and m~^ one can write down a 
representation for the effective Hamiltonian in the form 


Has=i — Hq + Hi. 


( 3 ) 


The hrst addendum Hq in the rhs of eq. (3) corresponds to the leading contributions in m~^ 
and coincides with the rhs of eq. (1) while the second addendum Hi is the leading order 
m~^ correction that comes from an annihilation diagram with the charmed quark running 
inside the loop. If one restricts the analysis to a scalar colorless gluon conhguration 
the additional contribution reads [12] 




t) - -UlsSfidr—G 

^ 120m2 " " ^ TT 


a 


( 8 ) 
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To derive an effective Lagrangian involving Goldstone bosons only, which can be used for 
the calculation of decay amplitudes, one has to replace a QCD operator 

= msSudiG^ (9) 

in eq. (8) with its mesonic counterpart. The following representation holds 

msSndiG^ = Af^{U^x )23 + Bf^{U^x) 23 tr{d^U^d^U) + other 0{p‘^) terms (10) 

where y stands for the meson mass matrix, U = jg the unitary matrix describing 

the octet of pseudoscalar mesons, A and B are dimensionless parameters to be computed. 
The second term on the rhs of eq. (10) is separated from other O(p^) structures for the 
reason that will be clarihed below. 

The O(p^) term in eq. (10) is a tadpole. The appearance of such a term is the 
consequence of working with a wrong solution for the ground state. This term merely 
renormalizes the effective Lagrangian of strong interactions and can be absorbed into the 
meson mass matrix by a suitable SU{3)l ® SU{3)r rotation [3,14]; it does not affect any 
observables. A contribution to the physical amplitude is determined by the 0(p"^) part of 
the chiral representation (10). The most transparent way to obtain this contribution is 
to consider the quark-gluon operator (9) as a product of the (pseudo)scalar quark current 
and the scalar colorless gluon operator. Then one can replace the quark operator by its 
mesonic realization according to the PGA G hypothesis [13] 

msSRdL -4^(G^x)23- (11) 

O 


On the other hand there is a low energy theorem based on fundamental properties of the 
energy-momentum tensor that gives the chiral representation of the gluon operator [15] 

+ o(p‘‘) ( 12 ) 

where 6 = 9 for three light quark flavors. Eqs. (11,12) give for the B parameter 



(13) 
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This approximation corresponds to the simplest physical picture where the kaon is annihi¬ 
lated by the pseudoscalar quark current while the pion pair is born by the gluon operator. 
Eqs. (8,10-13) lead to an effective chiral Lagrangian of the form 

_JLy A/* n _ 


it// = - T)y±::^yx)2Mayd'‘U). 


( 14 ) 


4806m^ 

Since the pion pair is born by the gluon operator this Lagrangian describes the decay 
channel with gluons forming an intermediate state. 

Now the corresponding K ^ nn decay amplitude becomes explicitly calculable. We 
use the standard parametrization of the amplitude Aq with the isospin transfer AJ = 1/2 


Gp 

ReAo ^ -^sin0cCosOcfKm^K9i/2 


( 15 ) 


where 6c stands for Cabibbo angle and gi /2 is a dimensionless parameter. The new 
contribution reads 


^9i/2 — 


m 


SObm"^ 


^ ~ 10“^ 


while experiment gives [16] 

9T/2 = 3 . 9 , 

and the most recent theoretical estimate is [2] 


(16) 


(17a) 


gi/2 ~ 2.6. 


(176) 


Thus the local (perturbative) part of new decay mode is negligible according to the general 
estimate of a scale of the leading order charmed quark mass corrections [12]. 

The local effective Hamiltonian (8) does not exhaust the whole physics of the meson- 
gluon transitions. It can not account for the long-distance contribution connected with 
the propagation of a soft u-quark round the loop of the annihilation diagram. Because 
of the lightness of the u-quark this contribution can not be represented as a local vertex 
and ultimately depends on the infrared properties of QCD. Its investigation requires some 
nonperturbative approach. 

In so doing we start with a tree level Hamiltonian which after decoupling of the c-quark 
has the form 

Hls.i = ( 18 ) 
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The quantity of interest is an effective theory realization of that part of the operator Q 2 
which is responsible for a transfer of the kaon into gluons. Invoking the results of our 
previous consideration we can write down this part in the following form 

= g^f^{U^x)23tr{d,U^d^U) (19) 

where g'^ is a dimensionless parameter. We should note that the chiral representation of 
the whole operator Q 2 contains a large number of structures but we are interested only in 
the part corresponding to the transition with the gluons forming an intermediate state. 
Thus the problem is reduced to computing the chiral coupling constant g^ that can be 
done by studying the appropriate Green’s function (GF) via QGD sum rules technique. 
For the technical reason working with a two point GF is preferable. In the given decay 
channel the pions are born by a gluon cloud therefore the gluon operator can play 

the role of an interpolating operator of the pion pair. Thus, it is natural to choose GF in 
the following form 

G(p) = l{0\T^G‘^Gl,(x)Q,(0)\K’’{q))e‘'^dx\,.„. (20) 

A remark about the chiral limit for the kaon in eq. (20) is necessary. The representation 
(19) hxes the correct O(p^) chiral behavior of the considered decay amplitude and does 
not depend explicitly on the kaon momentum. Keeping a non-vanishing kaon momentum 
leads to a shift of the decay amplitude that lies beyond the accuracy of the present 
approach. Thus we can put g = 0 in eq. (20) and work with GF depending on one 
argument only. 

Saturating GF in eq. (20) with 7r’''7r“ and 7r°7r° states (the lowest states with proper 
quantum numbers), substituting the Q 2 operator and using the low energy theorem (12) 
one obtains the following physical representation 

G’”'{p)=g<^^rfp*\n{-^) + 0{p‘). ( 21 ) 

TT^OJTT jjL^ 

The theoretical side reads after making use of OPE 

9 

X — T) (y. 

^ + 0(p“^). (22) 
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Factor rus in eq. (22) provides the correct chiral property of GF and justihes the represen¬ 
tation (19) for the operator Q 2 . By reducing the kaon state, one can transform eq. (22) 
into the following expression 


— (23) 

jj,^ TT 

where uiq determines a scale of nonlocality of the quark condensate, {qgsGf^i^a^^q) = 
m^{qq), ^^^(l GeV) = 0.8 ±0.2 GeV^ [17]. For extracting information about the chiral 

coupling constant we use hnite energy sum rules [18] with the result 

G ^ 3& flmlasiso) 

^ 128 TT ■ ^ ^ 

To take into account strong interactions at short distances the operator Qf in the 
effective Hamiltonian has to be multiplied by the corresponding Wilson coefficient Z 2 {sq). 
Finally, the new contribution to the theoretical estimate of the AJ = 1/2 amplitude in 
terms of the parameter gi /2 takes the form 


A5'i/2 


2^2(50) 


^^(so) 


TT 


8^2 


(25) 


This result needs some comments: 

1. This next-to-leading in l/iV/ expansion contribution is missed within the factorization 
framework and also within any approach where quark currents in four-quark operators 
are replaced by their mesonic counterparts separately. 

2. The gluon cloud in the intermediate state does not form a resonance state and, there¬ 
fore, the contribution (25) is not suppressed by a large scalar meson mass. 

3. In general, some more complicated scalar colorless gluon conhgurations, for example, 

appear as intermediate states in this channel as well. However the 
theorem (12) shows that the two pion form factor of such conhgurations could be of the 
O(p^) or higher order in chiral expansion that leads to a negligible 0(p®) shift of the decay 
amplitude. 

The question now is what numerical value for the duality interval Sq has to be used. 
Actually, the allowed value of the duality threshold is quite restricted by the form of the 
physical spectrum and by the requirement of absence of uncontrollable corrections. To 
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suppress contributions of higher mass states, for example, a scalar meson (t(0.9 GeV), to 
the considered channel one has to take Sq < (0.9 GeV)^. At the same time the physical 
representation (20) is obtained in the leading order in chiral perturbation theory and 
the whole procedure is justihed until the chiral expansion parameter So(87r^/^)“^ remains 
small [13,19]. On the other hand at the scale ju less than 0.8 GeV the perturbative 
corrections to Wilson coefficients become uncontrollable [2] and for consistency of the 
approach one has to set the low limit of the duality interval to be sq > (0.8 GeVY. The 
reasonable choice for the duality interval now reads sq = (0.8 GeV)^. 

Let us estimate the uncertainty of our result. On the physical side of sum rules the 
errors related to higher order terms in chiral expansion, which have been omitted in eq. 
(21), are, in general, unknown. But one can hope that in the spirit of chiral perturbation 
theory they are about 25% [13,19]. On the theoretical side of the sum rules the errors come 
from two sources. The first one is the perturbative part of OPE (the unit operator) that 
is suppressed by an extra loop factor ~ 10“^ and can not lead to a sizable change of 

our result. Next nonperturbative corrections due to operators with higher dimensionality 
seem to be more important. They start with the dimension eight operators which have 
already been discussed. Numerical estimates would require knowing the matrix elements 
of those operators between the kaon and the vacuum state which are not available now. 
But as a hrst approximation the relative weight of these corrections can be represented 
by the ratio 


) r(3) 1 

^ ~ 10“^ 26 

somg r(5) 

where ~ (330 MeV^ and the factor r(n) comes from a quark loop with 

n — 1 gluon held or mass insertions. The situation here is quite similar to that of the 
analysis of l/m^ corrections [12] where a contribution of dimension eight operators is 
suppressed rather numerically than parametrically. Taking into acconnt the nncertainty 
of determination of the parameter mg we estimate the error bonnd to be abont 40%. 
Numerically one obtains 


Agi/2 = 0.56 ±0.22 


( 27 ) 


at the point Sq = fV = (0.8 GeV^, ^qcd = 300 MeV, Z 2 {so) = 1.49 [9]. 



Thus, the new contribution provides about 15% of the observed amplitude (17a). At 
the same time it is comparable with the leading order result for the decay amplitude ob¬ 
tained by the naive factorization of the four-quark operator Q 2 when all strong interaction 
corrections are neglected 

dn = 5/9- (28) 

This implies the strong violation of factorization in the O(p^) order in chiral expansion 
that leads to additional enhancement of the theoretical estimate of the K tttt decay 
amplitude with the isospin transfer AJ = 1/2. 

To conclude, in the present paper we consider a new K tttt decay channel with 
the gluons forming an intermediate state. Both short-distance (perturbative) and long¬ 
distance (non-perturbative) parts of the corresponding amplitude have been calculated. 
The latter being dominant gives an additional enhancement of the K tttt decay with the 
isospin transfer AJ = 1/2 and provides about 15% of the of the experimentally observed 
amplitude. New contribution is of the 0{p^) order and is lost within the factorization 
framework. It allows us to conclude that there is a sizable violation of the factorization in 
the O(p^) order in chiral expansion. Taking into account corrections of this order along 
with the usually considered corrections to coefficient functions of Wilson expansion 
can help to resolve the AJ = 1/2 problem within the SM. 

This work was supported in part by Russian Fund for Fundamental Research under 
Contract No. 93-02-14428, by Soros Foundation and by Japan Society for the Promotion 
of Science (JSPS). 
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